August 2000 



Active-Sterile neutrino oscillations and BBN+CMBR constraints 



P. Di and R. Foot^ 

^ Istituto Nazionale di Fisica Nucleare (INFN) 
^ School of Physics 
Research Centre for High Energy Physics 
The University of Melbourne 
Victoria 3010 Australia 
( dibari,foot@physics.unimelb. edu. au ) 



Abstract 

We show how active-sterile neutrino oscillations in the early Universe can 
play an interesting role in explaining the current observations of CMBR 
anisotropics and light element abundances. We describe different possible 
phenomenological scenarios in the interpretation of present data and how 
active-sterile neutrino oscillations can provide a viable theoretical framework. 
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I. INTRODUCTION 



The standard big bang model is a simple and testable theory of the evolution of the 
Universe. One of the quantitative tests of standard big bang cosmology lies in its predictions 
of the primordial abundance of light elements. Standard Big Bang Nucleosynthesis (SBBN) 
contains essentially just one free parameter, the baryon to photon ratio at the time of 
BBN 7], while it predicts the primordial nuclear abundances of several light elements. Each 
measurement of a primordial abundance, that has to be inferred from observations at the 
present time, provides in principle a measurement of 77. The success of the theory relies on 
the consistency of the values for t] that can be deduced from different nuclear abundances. 

In recent years two different measurements of the Helium abundance have been given. 
A first group finds Q 'low' values Yp = 0.234 ±0.003, while a second group finds 'high' 
values Yp = 0.244 ± 0.002, which are mutually compatible at about 2.5 a level only. 

The SBBN provides numerically a relation Yp{rj). A linear expansion gives 0: 

Yp '''''' iv) = 0.2467 + 0.01 In . (1) 

In this way for the low value of Yp one obtains t^s^btv — 1-5 =t 0.4, while for the high value 
of Yp one obtains rj^^^j^ = 3.9 ± 0.8. 

Meanwhile, the value (D/H)^ = 3.39 ± 0.25 | has been deduced ^ for the primordial 
deuterium abundance from observations toward two high redshift quasars. Deuterium is 
an ideal 'baryometer' which gives an accurate measurement of baryon abundance in the 
context of SBBN. Indeed one finds {Vlb h'^)%%N = 0-019 ± 0.0024 (95% cl) § and from the 
simple relation rj ~ 273 ^Ib h'^, this corresponds to rj^sBN = 

5.2 ±0.65 (95% cl) []. This value 
is clearly not consistent with the r] value obtained from low Helium values although it is 
consistent with the rj value from the high Helium value. 

On a new front, two balloon experiments provided the first accurate measurements of 
acoustic peaks in the Cosmic Microwave Background Radiation (CMBR) anisotropics PJTO 



and from these observations it has been possible to infer a value for the baryon to photon ra- 
tio. The BOOMERanG experiment finds {VL^ h^)cBR = O.OSGto fot 10 ^^il^ the MAXIMA 



experiment finds {^bh^)cBR = 0.03llo;oo6 IQ- These two independent measurements are in 
quite good agreement and seem to exclude the presence of large systematic errors. A com- 
bined analysis of the two gives the result {flbh'^)cBR = 0.033 ± 0.005 [1^ that corresponds 



^Where not otherwise indicated, all errors are meant at 68% cl. 

neutron life time t„ = 887sec has been used. Here and everywhere r] is expressed in unit of 
■j^g-io_ rpj^g central value for r] = 5 has been updated according more recent analysis Note that 
this expression is accurate to within 0.001 for 3 < ry < 10 while is less accurate for r/ < 3. 

HD/H)5 = 10^ (D/H). 

In a very recent analysis even smaller errors are found: i^bh^)sLBN = 0-019 ± 0.0018 (95% cl) 
and correspondingly %^^jv ~ ^-^ ^ 0.5 (95% cl) 
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to rjcBR = 9-0 ± 1.4. This value is higher than the BBN predictions (given above) from 
the inferred values from both Deuterium and Helium. These discrepancies may be due to 
systematic errors but it is also interesting to consider possible explanations in terms of non 
standard physics. One possibility is that BBN and CMBR are probing different quantities, 



as they involve different physical mechanisms and at different times (see for example 
Instead we will consider this discrepancy as a hint for non standard BBN. 
We will consider two view points: 

1) The discrepancy is between CMBR and Helium while the discrepancy between CMBR and 
Deuterium is due to systematic uncertainties. This is plausible because it is very difficult to 
identify 'clean' absorption systems providing reliable measurements for deuterium and the 
quoted results that we used were derived from only two such measurements. 

2) The discrepancy between deuterium and CMBR as well as the discrepancy between 
Helium and CMBR are both real and due to non standard physics 0. 

The main purpose of this paper is to explore the possible explanations for these discrepancies 
in terms of active-sterile neutrino oscillations in the early Universe. 

At the present time there is very strong evidence for neutrino oscillations coming from 
atmospheric, solar and the LSND experiment (for a review, see e.g. Ref. [jrB|). The atmo- 



spheric neutrino anomaly can be solved (most simply) via approximately maximal z/^ Vr 
oscillations or via i/^ Ug oscillations (where represents a hypothetical sterile neu- 
trino). On the other hand, the observed solar flux deficit (about 50% of the expected value) 
suggests approximately maximal Ue u^^r oscillations or approximately maximal Ue — > i/g 



oscillations (see e.g. Ref. |T8| and references there-in). Finally, the LSND experiment im- 
plies the existence of small angle z/g oscillations with 5m? ~ leV"^. The combination of 
these three neutrino anomalies suggests the need for at least one sterile neutrino. Perhaps 
the most elegant solution to these neutrino anomalies poses that each neutrino is maximally 
mixed with a sterile partner (with small mixing between the generations) []T^. Of course, 
there are many other possibilities. In any case, for illustrative purposes we will focus on the 
simple case of z/g oscillations in isolation and discuss some of the other possibilities 

qualitatively where appropriate. 

Ordinary-sterile neutrino oscillations have remarkable implications for the early Universe. 
In particular, ordinary-sterile neutrino oscillations can generate large neutrino asymmetries 
in the early Universe pO|-p^ (see also Ref. [^]), so large in fact as to imply significant 



modifications to BBN In the simple scenario of oscillations in isolation, large 

Ve asymmetry is directly produced, while, in three (or more) meutrino mixing scenarios an 
a-neutrino asymmetry [a = n, r) could be first generated by z/q, Ug oscillations and then 



converted into an electron neutrino asymmetry by a z/q, — z/g oscillations |22,24]. Focusing on 



the simple case of direct production, we will find the values of mixing parameters which might 



^Note that we will neglect from our analysis measurements of ^Li abundance. This because in 
recent years various analysis conclude that there is still a big uncertainty on the level of depletion 
of the primordial abundance to the values that we currently observe (for a recent review see [15|). 
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explain the possible discrepancies in SBBN. Interestingly, it turns out that the suggested 
parameter space implies a Vf. with mass ~ leV which is close to the current experimental 
bound. Furthermore it is also consistent with the measurements of the LSND experiment 
and thus can be potentially tested in the near future at mini-Boone. 

The outline of this paper is as follows. In section II we briefly discuss the relation 
between VL^ inferred from X-ray measurements in galaxy clusters and the value suggested 
by CMBR. Interestingly both are consistent with Vtm = 1. In section III we show how 
the z/g — >■ Vg oscillation generated L,^^ can reconcile the high rjcsR with the BBN Yp results 
(case 1 above). In section IV we briefly examine ways in which the discrepancies between 
Deuterium, Helium and CMBR may both be reconciled (case 2 above) . In section V we will 
present a new possible phenomenological scenario, in which large scale inhomogeneities in 
the nuclear abundances are admitted. Also in this case active-sterile neutrino oscillations 
may provide a viable theoretical model. We conclude in section VI. 



II. COSMIC CONCORDANCE OR DISCORDANCE? 

In this section we make some comments on the previous indications of a 'baryon catastro- 
phe' from X-ray measurements in galaxy clusters and how it might be related to the inferred 
values of the baryon number from CMBR observations. Recall it is possible to estimate the 
baryon to total mass ratio in clusters of galaxies from X-ray measurements [^] obtaining 



(at la) that il^h^/ilj^ = (0.05 ± 0.01) pg] (this is confirmed also by measurements based 
on the SZ effect that give fib h/flm = (0.06 ± 0.006) |2^). If one imposes that flm = 1 and 
using a lower limit on h > 0.5 finds immediately that flfyh"^ > 0.035, much bigger than the 
upper bound that is deduced in a SBBN, both from Deuterium and Helium abundance (but 
consistent with CMBR as we will discuss in a moment). This 'baryon catastrophe' in SBBN 
was 'solved' by assuming that Qm can be much less than 1 which implies that we need to 
give up the inflationary paradigm {Qo = 1) or to admit the presence of a large cosmological 
constant A such that Qm + f^A = 1- In this way, using the SBBN value, previously given, for 
Qfyh'^ from Deuterium abundance, one can infer a value for = 0.45 ± 0.15. This picture 
has been supported by the discovery of an acceleration expansion from SNe la |2^, that 



also points to the existence of a large cosmological constant term, ||2^- Roughly these 
measurements provide the constraint, in the flm — plane, 

^]A = 1.3^]^ + 0.4±0.2. (2) 

After the SN results, the first accurate CMBR measurements of the first acoustic peak 
position seem to confirm the idea of a flat Universe |3^. In this case one immediately 



deduces from Eq.(g) a value flm = 0.25 ± 0.1 in very good agreement with galaxy cluster 
measurements when the SBBN value for flfy h"^ is assumed: three independent methods match 
each other, with good accuracy, in a region around the point (0.3, 0.7), in the plane {flm, ^a) 
('cosmic concordance' PT[]). 

It is clear however that if one now uses the new CMBR estimation for fl^ h?, then from 
galaxy clusters one obtains that flm = 0.8 ±0.3 suggesting 'cosmic discordance' (albeit only 
mildly) between galaxy cluster measurements and SNe type la. Moreover now the existence 
of a cosmological constant is not required any more from galaxy clusters . This seems to 
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be suggested also from a recent analysis of CMBR data when a large neutrino asymmetry is 
allowed |33,32]. Thus, overall things are not so clear at the moment. It has also been argued 



in Ref. |Q that the Supernovae evidence that the expansion of the Universe is accelerating 
is not yet compelling. Whether or not a large cosmological constant exists needs to be 
confirmed independently, perhaps by future analysis from the planned satellite experiments, 
MAP and PLANCK will help. In the meantime, for people with the theoretical prejudice 
that is negligible, there is now some good news since X-ray measurements from galaxy 
clusters and the CMBR results are both consistent with = 1. 



III. CMBR ANISOTROPIES AND HELIUM OBSERVATIONS 

Even assuming the high value for the Helium abundance, the resulting value of rj and 
that one deduced from CMBR differ at about 3a level. 

Let us consider the situation from a formal point of view that will be useful for further 
developments. In the SBBN picture the experimental constraint ¥^^^^{1]) = Yp^^ gives 
a measurement of 77 = tfs^^BN^ observed in the introduction. However CMBR gives 

an independent measurement of 7] and, assuming it to be a reliable one, provides a simple 
test for SBBN as now one has to satisfy the constraint Y^^^^ {ticbr) = Y^^^. Present 
measurements do not pass this test and thus SBBN is somewhat discrepant with Helium 
and CMBR observations. 

Taking this as a hint for new physics, it suggests that we need to modify SBBN intro- 
ducing a new parameter X. It is clear that, if in the physical ranges of values for X, all 
values of Yp are possible (with r] = ticbr), then it is always possible to find a value for X 
satisfying the test Yp^rjcMBR, X) = Y^^^ in the non standard BBN model. 

It has been known for a long time that this 'game' can be performed allowing a modi- 
fication of the standard particle content before the BBN epoch |3^. This modification can 



be parametrized with the (extra) number of (light) neutrino species [BB 



AiV.^ = E - 3, with N^^^/-^, (3) 



{X = z/g, z/^, z/t- + new particle species), where Ti, = TdRd/R is a fiducial temperature of 
ideal neutrinos that would instantaneously decouple at Td ^ me/2 (but also -C m^/2), 
without sharing any entropy release, from electron-positron annihilations, with photons 0. 

With this extra parameter, the SBBN prediction for Yp is modified and approximately 
the change is given by: 

AYpi,, AN^) ^ I r/---(,) (1 - Yl^^^ !Ihl_pL ^ ^ 0.012 AN^, (4) 



^Note that with this definition, in the standard model of particle physics one finds that ANp is not 
exactly zero, due to the fact that actually neutrinos are slightly reheated during electron-positron 
annihilations (see |^7j and references therein). 
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where T"'^^ ^ 0.75MeV is the freezing temperature of neutron to proton ratio in the standard 
case, while N^^^ = 43/8 is the number of particle species and again the last expression has 
been evaluated for rjcBR — 9. 

It has also been known for a long time that the standard prediction for the neutron to 



proton ratio is modified if one allows the electron neutrino and anti-neutrino distributions 
in momentum space to deviate from the standard case in which the thermal equilibrium 
distributions with zero chemical potentials are assumed Q. In this case an infinite number 
of new non standard parameters, the values of the extra numbers of electron neutrinos 
and antineutrinos in each quantum state for any momentum, can be virtually introduced. 
However, in realistic models, usually the deviations depend on a finite number of parameters. 

A particularly simple model is obtained when the distributions depend only on the 
neutrino degeneracy C,e, playing the role of a new parameter X. This has to be generated 
earlier than the BBN epoch and it is usually assumed that it is generated also before electron 
neutrino chemical decoupling so that + = 0. This corresponds to having a neutrino 
asymmetry given by: 



In this case the modification of the SBBN prediction [see eq.(|^)] for ^ 1 is given by: 

AY,ir^,U ^ -Yf^^^'iv) [l - ' ^ ^^' j ^ -0.22ee, (6) 

where the last expression has been calculated for ticmbr — 9. Thus values ^ 0.035 and 
~ 0.08 can easily solve the discrepancy between CMBR observations and high and low 
Helium values respectively. 

More generally one can allow a deviation from the standard prediction for the Helium 
abundance due both to a modification of the expansion rate and to distortions of the electron 
neutrino and anti-neutrino distributions. In this case one can distinguish two different 
contributions to the variation of the Helium abundance compared to the standard case: 

Arp = AF/ + AF/-. (7) 

Note that this distinction is not ambiguous as it could appear. One can in fact always 
calculate at any instant the quantity /S.NP from the eq. and from that deduce the 
corresponding value of defined as the value of AYp when the standard electron neutrino 
and antineutrino distributions are assumed. Afterwards one can calculate AY^'^^ = AYp — 
AYP. It will prove to be convenient to define also a total effective number of neutrino species 
ANy that combines both the effect of a modification of the expansion rate and that one due 
to the distortions of electron neutrino and anti-neutrino distributions: 



^Electron neutrino and anti-neutrinos distributions are directly involved in determining the rates 
of the /3-reactions (n + <-> p + e~ , n + e"*" <->■ p + S'e) responsible for the final value of the neutron 
to proton ratio together with the neutron decay. 
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AN, = ANP + AAT/-, with AAT/- = (8) 

With these definitions, the procedure to calculate AN^, and the specific contribution AN^"'' 
is particularly simple if one uses the linear expansion (^: AA^j^ ~ Ay^/0.012 and thus 
ANf-- ~ Ayp/0.012 - AA^^ Q. 

In the case of an electron neutrino asymmetry created before the electron neutrino chem- 
ical decoupling ^e + Ce = 0) it is easy to see that AN^"^ ~ — 18 |- In order to reconcile 
the discrepancy between the Helium abundance measurements and the CMBR observations, 
a negative AA^^ is required. More precisely, imposing the constraint Yp{rjcBR, ^Ny) = Y^^^, 
one finds: 

AA^^ = -1.5 ± 0.4 for F/^p = 0.234 ± 0.003, 

AA^^ = -0.7 ± 0.3 for F^^^p = 0.244 ± 0.002. (9) 

One could simply hypothesis that a large pre-existing asymmetry exists, however a more 
subtle (and testable as we will show) possibility is that light sterile neutrinos exist. Naively, 
such neutrinos might be expected to lead to a positive AA^j, as they would simply generate 



a positive AN^, however it has been shown ||2^J40|J24|| that due to the dynamical generation 



of by the oscillations themselves, there will be also an important contribution to AA"/''^ 
and the total effective AA"i, can be negative if Ly^ > f^. 

This contribution cannot be expressed in terms, for example, of the final asymmetry by 
a simple relation as in the case of a pre-existing asymmetry. This is because the asymmetry 
is generated, in the interesting cases, below the chemical decoupling temperature and even 
below the thermal decoupling temperature and thus is changing during the time near the 
freezing of the neutron to proton ratio. Moreover, as the thermal equilibrium assumption is 
not satisfied anymore, the electron neutrino distribution will deviate from equilibrium and 
this effect has also to be included. Thus the results can be only calculated numerically. 

The sign of L^^ cannot be predicted because it depends on the sign and magnitude of 
the initial lepton number asymmetries. For the purposes of this paper, we assume that it is 



®0f course the choice to describe results in terms of AYp or of AA'^^ is simply a matter of taste. 
One can say that in the first case the astrophysical point of view is more emphasized than that one 
of particle physics or vice versa in the second case. In this paper we clearly prefer the second one. 

^We remind that this is valid only for 1. Note also that for large neutrino asymmetries 

(^ej^^iCr ~ 0.5) these would also give a non negligible contribution to AN^ 

^'^The problem to get a number of effective neutrinos less than three is not new. It also arises in or- 
der to alleviate the tension between Deuterium measurements and low values of Helium abundance 
('BBN crisis') [^. Also in that case the same non standard solutions can be invoked. Rounding 
up the usual suspects, we have a MeVr neutrino decaying prior to the onset of BBN the 
existence of a large electron neutrino asymmetry |43| and active-sterile neutrino oscillations |22|, 
that we are re-considering in this new context. 
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positive since we need to generate AA^^, < 0. The simplest example of neutrino oscillation 
generated L^^ is the direct production of L^^ by z/g — >■ oscillations. In this case we can 
ignore the oscillations involving z/^, provided that either their masses are very small (so 
that the largest \5m?\ belongs to the z/g oscillations and the other oscillations have 

|5m^| much less than 1 eV"^) or that they do not mix with the z/g, (i-e. the z/g, z^s decouple 
from the z/^, z^r in the neutrino mass matrix). In this way the mixing is simply described 
by two parameters, the difference of squared eigenstate masses and the mixing angle in 
vacuum, sin^ 26*0. 

In Figure 1 we solve the quantum kinetic equations for z/g — > z/^ oscillations for sin^ 29q = 
10^^ and dm? /eV'^ = —0.25, —0.5, —1.0, —2.0, —4.0. (For details of the numerical procedure 



see Ref. [^). Let us now discuss the behaviour exhibited in this figure. As already discussed 
in detail in previous publications [p0|-p2| the evolution of lepton number can be separated 
into three distinct phases. At high temperatures the oscillations are damped and evolve so 
that L^^^ <^ 7] (where L^'^^ = 2Ly^ + Ly^ + L^^ +r], and r] is related to the baryon asymmetry). 
In this region the resonance momentum for neutrino oscillations is approximately the same as 
anti-neutrino oscillations. If < (which means that the mass eigenstate which is mainly 
Ug is lighter than the mass eigenstate which is mainly z/g) then at a certain temperature, Tg, 



which is given roughly by pU[ , 



T.-^^{=^^^) (10) 

exponential growth of neutrino asymmetry occurs (which typically generates a neutrino 
asymmetry of order 10~^, as shown in figure 1). Taking for definiteness that the Lj,^ is 
positive, the anti-neutrino oscillation resonance moves to very low values of p/T ~ 0.3 while 
the neutrino oscillation resonance moves to high values p/T ^ 10 (see Ref. [^] for a figure 
illustrating this). The subsequent evolution of neutrino asymmetries, which is dominated by 
adiabatic MSW transitions of the antineutrinos, follows an approximate 1/T^ behaviour until 
the resonance has passed through the entire distribution. The final asymmetry generated 
is typically in the range 0.23 ~ L^y^ ~ 0.37 Because the oscillations are dominated by 



adiabatic MSW behaviour it is possible to use a relatively simple and accurate formalism to 
describe the evolution of the system at 'low temperatures', T ~ Tc/2. In fact, we only need 
to know the values of the oscillation resonance momentum at T ~ Tc/2. Previous numerical 
work has already shown P2| that by T ~ Tc/2, neutrino asymmetry is generated such that 



0.2 ~ p/T ~ 0.8 (the precise value depends on sin^ OoySm^). Furthermore the subsequent 
evolution is approximately insensitive to the initial value of p/T in this range. 



For full details of the evolution of T^,^ and ANi^ in this model see Ref. . The evolution 
of the momentum distribution of electron neutrinos is also computed and fed into a BBN code 
(that is solved concurrently) which allows us to compute Yp for each choice of and sin^ 29o 



Particularly simple results are obtained when the constraint sin^ 26'o \/(|5m^|/eV^^ ^ 



2.5 X 10~^ is imposed. This corresponds to having AA"^ ~ 0.1 prior the onset of the neutrino 
asymmetry generation |^4|. Moreover, for the interesting values \5m?\ -C 100 eV^, most of 
the generated neutrino asymmetry and its associated sterile neutrino production, occurs 
below chemical decoupling so that AA^^ remains negligible. In this way the only significant 
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contribution to derives from the AN^"^ part which arises from the depletion of the 

i>e states as the MSW resonance passes creating Lj,^ in the process. For ~ —1 eV^, 
the large neutrino asymmetry is generated provided that sin^ 26q ~ few x 10~^° pl|j22 



(which is essentially the adiabatic condition for this system). With these two constraints 
on mixing parameters, the resulting AA^^^ is practically independent of sin 26o and thus we 
have a full correspondence AN^ ^ 5m^. The result is given in Figure 2. From this figure 
we can translate the constraint Yp{rjcBR, ANy) = Y^"^^ on AA^,^ [see Eq.Q, into a sort of 
'measurement' of (5m^, i.e. 



ANy = -1.5 ± 0.4, Sm^ = -2.5 ± 1.0 eV\ 

ANy = -0.7 ± 0.3, Sm^ = -0.8 ± 0.5 eV^. (11) 

These values of Sm"^ are interesting from several points of view. They imply that m^^ ~ 1 eV 
(assuming m^^ <^ n^'Ue) which is close to the present experimental limit. Furthermore, if 
is heavier than the z/^ state then the LSND , Smf^^^ is, approximately, the same as 
the (5m^ for Ue —>■ Ug oscillations. Thus, if this simple scenario is the cause of the BBN 
discrepancy it can be potentially tested in the near future at mini-Boone. 

Above we have discussed things in the model where L^^ is produced directly. It is also 
possible to produce L^^ indirectly. E.g. if Ur is the heaviest neutrino and oscillations between 
Ur — ^ generate a large L^^ some of which is transferred to L^^ by z/g — >■ Ut oscillations 
p^ , |40|p^ . The indirect mechanism typically generates a smaller L^^ leading to A A",^ in the 



range —0.7 ~ AA^,^ ~ if L^^ > 0. Models with three sterile neutrinos (such as models with 



mirror neutrinos) have also been studied ||4^. These models can also accommodate negative 
AA^j^ in the range —1.5 ~ AA^j, ~ if L^^ > 0. In fact it is fair to say that a deviation of 
AA^j, from zero is a generic consequence of models with light sterile neutrinos if one of the 
active neutrinos has mass in the eV range. 



IV. 'JUST SO' BBN ? 

We want now to include the Deuterium observations in our analysis. In this case the 
discrepancy between CMBR and nuclear abundances observations, becomes even more puz- 
zling. In the SBBN {D/H){r]) oc ?7~^-^ |^6| and this means that having tjcbr ~ 9 corresponds 



to (D/H)^ ~ 1.5, a quantity about half the measured one. One could hope that, within the 
model of z/g ^ oscillations with AA^^ <^ 1 discussed in the previous section, choosing the 
values of able to reconcile CMBR and Helium observations, it would also be possible 
to satisfy the constraint [D / H){ricBRi5Tn^) = {D/HY^p. This is however not the case as 
the negative values of AN^"^ leaves almost unchanged the standard value corresponding to 
VcBR- Thus the only way out is to enlarge the space of parameters in the model of BBN. 
This can be done allowing also a non zero AN^. This possibility has also been studied for a 
long time and recently reproposed in to solve the BBN-CMBR discrepancy. In this 



way the problem is now to find values of AA^^ and AA'/''^ that satisfy simultaneously the 
constraints: 

Y^ivcMBR. AN^, ANt) = F-P, (12) 
{D/H){'ncMBR. ANP, ANt) = (D/Hr^ (13) 
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In a recent analysis [^], in which a pre-existing electron neutrino asymmetry is assumed, the 
authors find that to reduce the discrepancy within a 2a level, a range of values 1 ~ AA^^ ~ 11 
and correspondingly 0.07 ~ 0.43 must be chosen. 

As we said in the previous section, it is not possible to make a straightforward comparison 
with the active-sterile neutrino oscillations, as the effect of the generation of a neutrino 
asymmetry is not easily related. However we can make some qualitative comments. It is 
quite easy to have marginal consistency at the 2a level by just modifying the constraints 
imposed on the mixing parameters to solve the discrepancy of CMBR with the Helium 
abundance alone. In fact simply increasing the mixing angle, with a fixed 5m?, such that 



sin2 26'o y (|(5m2|/eV^) ~ 2 x 10-^ one gets a AA^^ ~ 0.6. It is likely that the values for 
found in the previous section will be slightly increased as a higher AiV^''<= is now required to 
satisfy also the constraint from the Helium abundance |^. Therefore, within the framework 
of active-sterile neutrino oscillations, the search for the suitable values for AN^'''^ , AA^'' is 
translated in a search for the right 6m?, sin^ 29o values. 

It is clear however that allowing for the existence of just one sterile neutrino species, val- 
ues of AA^^ ~ 5, required to have a best fit, are not possible. In this case one has necessarily 
to assume the existence of more than one light sterile neutrino species. One amusing possi- 
bility is the idea that a mirror world exists where every particle has a corresponding mirror 
particle (see also Ref. [19| and references there-in). The main theoretical motivation 



for this theory is that it allows parity and time reversal to be exact unbroken symmetries 
of nature. In the context of this theory, it is usually assumed that the temperature of the 
mirror particles is less than the ordinary ones in the early Universe P3|,IT3. However, it is 



possible that interactions between the ordinary and mirror worlds may be strong enough to 
thermalize the mirror pRrticlcs sucti tti3jt TfdiyyQ'p — TQ^^difiary ^ • this case AA^^ ~ 6.14. 
To reconcile BBN with such a large value of AA"^ requires a large ~ 0.4 pre-existing 
asymmetry. (It needs to be pre-existing because if Tmirror = T, then there are equal den- 
sities of ordinary and mirror neutrinos which means that one cannot generate significant 
asymmetries). Alternatively, if 0.7T < Tmirror < T, then 1.5 ~ AA^^ ~ 6.14. In this case 
neutrino oscillations can generate a significant asymmetry which may potentially lead to 
a model consistent with BBN for a range of parameters^. 



^^It must be said that increasing the mixing angle there is a region where at the onset of the 
asymmetry generation rapid oscillations are found (see also Ref. [50|). It is still an issue 



whether this is a real feature if the solutions or just simply an effect due to numerical inaccuracy. 
However if this effect would really exist, it is possible that the sign of the asymmetry could be 



randomly determined in different points of the space with the creation of lepton domains |51|. 
This would spoil the effect that we want to get, as in this case negative values of AA^^ would not 
be allowed. In any case at the high values of mixing angles that we are requiring in order to have 
AA^^ ~ 1 (of course AA'^^ cannot be too close to one otherwise this will suppress the final neutrino 
asymmetry), there are surely no rapid oscillations p9| . 

It may also be possible for Tmirror > T, with large L^^ generated by u'^ — > oscillations, 
leading to a consistent model for a range of parameters. 
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While such possibihties are interesting and testable, it may however seem surprising that 
nature should have large AN(^ and AiV/''<= which roughly cancel 0. Maybe the discrepancy 
will be alleviated from more precise measurements of t] from CMBR and a mild compensation 
with one or two extra neutrino species and a not too big neutrino asymmetry would be 
perhaps reasonable. From this point of view a crucial test in the future will be provided 
when CMBR will also be able to measure Nf^, while at the moment it only provides a rather 
poor upper limit iV^ < 13 |B3. It is however possible to imagine a different kind of solution 



within non standard BEN models that circumvent the requirement of a fine tuned solution. 
We now turn our attention to one idea in this direction. 



V. INHOMOGENEOUS NUCLEAR ABUNDANCES ? 

CMBR measures the baryon abundance on the whole observable universe with a co- 
moving size of GOOOMpc/i"^ and this would correspond to a SBBN prediction of 1^ ~ 0.25 
and {D/H)^ ~ 1.5. The Deuterium measurement is deduced from two Lyman absorption 
systems at Zabs — 3, corresponding to comoving distances of about 2000Mpc. The size of 
these systems is approximately equal to the comoving size of galaxies (lOOKpc — IMpc). 
Primordial Helium abundance values are deduced from ionized gas surrounding hot young 
stars at distances within ~ lOOMpc around us. It is then possible to imagine that an inho- 
mogeneous electron neutrino asymmetry could be the reason for the apparent discrepancies 
between Deuterium and CMBR as well as between Helium and CMBR. The quantities 
Yp^^^ {rjcBR) — 0.25, D / H^^^^ {ticbr) — 1-5 provide us the values of the nuclear abun- 
dances as they would be in absence of neutrino asymmetry. This would mean that in the 
absorption systems that we observe, a large negative neutrino asymmetry is needed to change 
(D/H)^ from ~ 1.5 to ~ 3. On the other hand to explain values of Yp less than 0.25 in our 
surroundings, as we already discussed at length, a positive neutrino asymmetry is required. 
Note that a hint of the presence of inhomogeneities in Deuterium abundances comes from 
the observation of Deuterium in a 2;abs = 0.701 toward QSO 1718+4807 where it was found 
{D / H)^ = 25±5 [^]. Other authors repeated the analysis and, even though they confirm an 



high value, they arrive at a much looser bound, {D / H)^ = 8 — 57, concluding that the deter- 
mination of D/H from QSO 1718+4807 is uncertain |^ |^. Using a more elaborate model 
for the velocity distribution inside the absorber, a third group finds (D/H)^ = 4.1 —4.7 []57| , 



in any case still higher than the value, 3.3 ± 0.25 deduced from the two cleanest absorption 
systems. Another system gives a result {D/H)^ < 6.7 |ES|. In a recent review the possibility 



of high amplitude inhomogeneities with an equal proportion of low values [D / H)^ ~ 3 and 
high values {D / H)^ ~ 10 is excluded ||T^. However it cannot be excluded that rare peaks 



model employing a decaying Mev r neutrino has also been recently proposed to get a 'just 
so' BBN scenario 

^^They observe in fact that in this case the spectra of the Lyman series lines is missing. This is 
needed to determine the velocity distribution of the Hydrogen and these measurements with the 
high value assume a single velocity component. 
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with {D/H)^ ~ 10 are present and in any case inhomogeneities with values changing in the 
range (D/H)^ = 1 — 4 cannot be excluded at the moment. The possibility for Deuterium 
abundance inhomogeneities has already been explained with the presence of inhomogeneous 
electron chemical potential an interpretation that could be now enforced by CMBR 



data. 

Active-sterile neutrino oscillations can give rise to an inhomogeneous field of electron 
neutrino asymmetry when the presence of small inhomogeneities in the baryon number is 
assumed In this case, the generated neutrino asymmetry can have an inverted sign 



in points where the baryon number is lower than the average value. Large scale inhomo- 
geneities in the electron neutrino asymmetry might be expected to generate inhomogeneities 
in the energy densities that would leave an imprint in the CMBR anisotropics that we do 
not observe. However in the case of active-sterile neutrino oscillations, inhomogeneities in 
the electron neutrino asymmetry would be compensated by inhomogeneities in the sterile 
neutrino asymmetry in a way that the energy density remains homogeneous and the mech- 
anism is not constrained by CMBR. There is one difficulty however due to the fact that 
one has to require the simultaneous presence of large scale regions with positive electron 
neutrino asymmetry and negative neutrino asymmetry. It has been shown that do- 
mains with inverted sign bigger than lOKpc cannot be generated. In this case even though 
at the onset of BBN one would get values of abundances in regions with positive electron 
neutrino asymmetry and also in regions with negative neutrino asymmetry, later on astro- 
physical processes, like supernovae explosions, would mix the different elements leading to 
approximately homogeneous values for the abundances. 

One way to circumvent this is to assume the existence of two scales. On small scales 
(less than the diffusion length at the time of freezing of neutron to proton ratio, ~ lOOpc) 
baryon number inhomogenities have to be present. On large scales, as big as required by 
Deuterium observation, the amplitude of these inhomogeneities has to change such that only 
in the regions where it is large enough a structure of small scale lepton domains with both 
signs can form. Neutrino diffusion would afterwards make them merge such that, in these 
regions, electron neutrino asymmetry is diluted to negligible values prior to the freezing of 
the neutron to proton ratio. On the contrary in the regions where domains did not form, a 
non zero electron neutrino asymmetry, with the normal sign, would be present []60 |. 

In this way one can easily get a field of neutrino asymmetries with values changing 
between zero and some maximum values. This means that it would be easily possible 
to accommodate CMBR with only Deuterium observations (in this case the normal sign 
should be positive) or with only Helium observations (in this case the normal sign should 
be negative). 

If we want to accommodate both Deuterium, Helium and CMBR then we need domains 
with both positive and negative signs on scales larger than 10 Kpc. However, as we men- 
tioned earlier, this violates the bound from Ref. [Q. Actually, the conclusion that domains 
with inverted sign on scales larger than lOKpc cannot be obtained relies on a simplified as- 
sumption for which domains with inverted sign cannot merge with each other. This is what 
would happen in the presence of a simple spectrum of baryon inhomogeneities with just two 
characteristic scale lengths as we just described : one is the scale of small lepton domains 
(^ lOOpc) and one is a scale that modulates the amplitude of baryon inhomogeneities in a 
way that in some regions lepton domains can form and in some others cannot. However in 
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a realistic spectrum of baryon fluctuations with a random presence of Fourier components 
things can be much different and one cannot exclude a priori that in some regions, small 
inverted sign lepton domains can occupy most of the space and they can merge with each 
other to form very large domains with a scale higher than about lOOKpc, both with positive 
and negative sign. These very large scale domains would give rise to inhomogeneities in 
the nuclear abundances that could not be washed out by astrophysical processes and would 
survive until the present. Q A clear signature of this mechanism would be the detection of 
high values of Helium (~ 0.30) in the regions where Deuterium is also measured with a value 
[D / H)^ ~ 3 while if some peaks with {D/H)^ ~ 10 really exist, here the Helium abundance 
should be even at level of 1^ ~ 0.50 [Q. However these measurements at large distances 
seem, at the present, to be quite challenging. Anyway when more measurements from Lyman 
absorption systems will be available, a clear signature of inhomogeneities could be possible. 
On the other hand if observations will exclude Deuterium abundance inhomogenities in the 
range (D/H)^ = 1 — 4 or larger, then an explanation of the BBN-CMBR discrepancy in 
terms of a spatially variating electron neutrino asymmetry, as we are proposing, would be 
ruled out. 



VI. CONCLUSIONS 



We have discussed the discrepancy between the inferred baryon number density, rj from 
recent CMBR measurements and the value inferred from standard big bang nucleosynthesis. 
This discrepancy may be due to some type of systematic error or may hint at new physics. 
We have explored one possible explanation in terms of active - sterile neutrino oscillations. 
We focussed on the simplest example to illustrate this possibility, and that is the direct 
production of L^^ by Ue Vs oscillations. Within the context of this model, we have 
shown that 5m? —\ eV"^ is required to solve the discrepancy between CMBR and Helium 
measurements and this would suggest that the electron neutrino mass is about 1 eV. This is 
a particularly interesting value, since it is right near the boundary of current experimental 
measurements. While we focussed on the largest discrepancy between r/^^SAf rjcBR, 
we also discussed the %BBAf tjcbr discrepancy, and its implications for models with 
sterile neutrinos. In particular, we looked at two possible scenarios. The first one, would 
reconcile the deuterium discrepancy with a large AiV^, while still needing a AN^^'^ of the 
opposite sign to reconcile the Helium measurements. We also proposed a second scenario in 
which we argue that an inhomogeneous electron neutrino asymmetry could exist which solves 
these discrepancies. For both of them we showed how active-sterile neutrino oscillations can 
provide a viable theoretical framework. 



^^We have to mention that in another alternative model proposed in [59| the simultaneous pres- 
ence of regions with positive neutrino asymmetry together with regions with negative neutrino 
asymmetry is a natural consequence. Here of course we are concentrating our attention on active- 
sterile neutrino oscillations, but the consideration that CMBR could be pointing to the presence 
of large scale inhomogeneities in the nuclear abundances has a general validity. 
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Clearly things will soon become more interesting as more accurate measurements of 
CMBR and light element abundances are done, and also, as we learn more about neutrinos 
from current and future experiments. Thus, it seems that large neutrino asymmetries, as 
generated from active - sterile neutrino oscillations offer an exciting interconnection between 
the rapidly developing fields of neutrino physics and early Universe cosmology. 
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Figure Captions 

Figure 1: Evolution of L^^^ = 2L,y^ + ri for z/g Us oscillations with sin^ 20q = 10~® and, from 
left to right, Sm'^/eV'^ — —0.25, —0.5, —1.0, —2.0, —4.0 obtained from numerically solving the 
quantum kinetic equations. The initial L,^^ = is taken and rj — 5 x 10~^° is assumed. Of 
course the low temperature evolution is approximately independent of these values. 

Figure 2: Change in the effective number of neutrinos for BBN, AiVj, versus —drn^ for the 
case L^^ > 0. 
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